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1. Introduction 



As is well known in standard plane-wave impulse approximation (PWIA) the differ- 
ential cross section for A(e, e'p)B reactions factorizes into an elementary cross section 
a ep , describing electron proton scattering, and a spectral function S(E m ,p m ), describing 
the probability for finding a proton in the target nucleus with energy (E m ) and momen- 
tum (p m ) compatible with the kinematics of the process. The differential cross section in 
PWIA is thus written as 

d 6 a 



de' ' dVt e dE p dVt p 



X<J ep S(E m} p m ) (1) 



with x a kinematical factor. It is this factorization property that makes A(e, e'p)B 
reactions so appealing for investigations of nuclear structure. 

According to eq. ([!]) single-particle distributions can in principle be probed in great 
detail. In fact several orbitals have been mapped out in many nuclei over rather extended 
momentum ranges in spite of the limitations of the PWIA approximation. These 

limitations are also well known. Final-state interactions and, more generally, distortion of 
both electron and nucleon wave functions due to electromagnetic and strong interactions 
with target and residual nuclei destroy the elegant simplicity of eq. ([[]). 

Nevertheless this equation is still very useful and is the basis for interpretation of 
experimental data. The latter are usually compared to theory defining an effective spectral 
function 

~ ! d 6 a 

S (E m ,p m ) = ( X a ep ) — ^ ^ (2) 
and an effective proton momentum distribution (or reduced cross section P, [§) 



P iPm) = S (E m ,p m ) dE m , (3) 

JAE m 

that in PWIA corresponds exactly to the density in momentum space of the single-particle 
orbital associated with a selected value of the missing energy E m . In general S and P 
depend not only on p m and E m but also on the remaining kinematical variables. 

Comparisons of theories taking into account in various ways final-state interactions 
and/or Coulomb distortions are also made for this effective momentum distribution. When 
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final-state interactions and Coulomb distortions are fully taken into account, departures 
between experiment and theory can solely be attributed to limitations of the single- 
particle nuclear model and impulse approximation. Thus eqs. ([!]) to (0) form also a basis 
for the empirical study of nucleon-nucleon correlations in the target. 

In recent years much work has been done along these lines on both experimental and 
theoretical fronts. Experimental exclusive (e, e'p) measurements have been made with 
high precision on quite extensive missing momentum regions @ -[[§]. The analyses of data 
based on standard distorted-wave impulse approximation (DWIA) have met two major 
difficulties. On the one hand the spectroscopic factors extracted from DWIA analyses 
of low p m (p m < 300 MeV) data are too small compared with theoretical predictions. 
As an example the extracted occupations of 3si/2 and 2d§/2 orbits in 208 Pb are S a ~ 
0.5, while theories on short-range correlations || predict at most a 30 % reduction of 
mean field occupations for levels just below the Fermi level. On the other hand the 
high-p m data (300 MeV < p m < 600 MeV) on p (p m ) for the same levels are much 
larger than the results of DWIA calculations that are compatible with those that fit 
the low-p m data Q. Although short-range correlations are expected to increase the 
high momentum components, their effect is negligible || at the small missing energies of 
the existing high-p m data, and effects of long-range correlations have been invoked 
The above mentioned DWIA analyses are based on non-relativistic independent particle 
models and use phenomenological (real and imaginary) Woods-Saxon potentials with 
parameters fitted to individual nuclei and shells. 

Alternatively, in recent years, the relativistic mean field approximation has been suc- 



cessfully used for the analyses of both low- and high-p m data [|7fl-fl0|l. In the relativistic 
distorted-wave impulse approximation (RDWIA) bound and scattered nucleons are de- 
scribed by solutions of the Dirac equation with scalar and vector (S-V) potentials. In the 
work of refs. |8[|-[10[| the scattering nucleon wave functions were obtained solving the Dirac 



equation with the relativistic optical potentials of ref. |TTJ and the bound nucleon wave 
functions were obtained with the TIMORA code ]T2[ based on Walecka's mean field ap- 



proximation [y| . The parameters of the S-V potentials are obtained from global fits and, 
at variance with DWIA, in RDWIA there were no parameters fitted to the particular levels 
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other than the spectroscopic factors. In all the cases studied the extracted spectroscopic 



factors are larger than in DWIA and are valid for low- and high-p m data ||7|-|10fl. For 
instance, for the above mentioned 3si/ 2 and 2d 3 / 2 shells in 208 Pb values of S a ~ 0.7 have 
been obtained J7L ffl , consistent with theoretical predictions, and a reasonable agreement 



between RDWIA calculations and experiment has been found in the high-p m region [10]. 

A global comparison of the unsealed results shows that, compared to DWIA, the 
reduced cross sections obtained in RDWIA are smaller at low p m and are larger in the 
high-p m region, thus providing a consistent description of all available data with moderate 
values of the spectroscopic factors. 

The reasons why RDWIA produce smaller cross sections in the low-p m region were 
investigated in refs. HQ. In ref. |J the success of the relativistic analyses was traced 
back to the improved treatment of distortion effects in the electron and outgoing proton 
waves. In the high p m region other factors can be important. In particular, the presence 
of higher momentum components in the relativistic bound nucleon wave function, and in 
the relativistic nucleon current operator can be expected to play a role in this region. In 
view of the success of RDWIA calculations it is worthwhile investigating more deeply the 
dependence of the results on the lower components of the bound nucleon wave function 
and on the choice of the current operator. This is the purpose of this paper where we 
restrict ourselves to the relativistic plane-wave impulse approximation (RPWIA). 

The reasons why we consider here only the plane-wave limit are twofold. On one 
hand it is convenient to restrict the study to RPWIA to be able to disentangle effects of 
distortion from effects of high p m components in the bound nucleon wave function and 
current operator. On the other hand it is interesting to study this approximation be- 
cause an important difference between relativistic and non-relativistic approaches already 
appears in the plane- wave limit. Namely, the factorized expression in eq. (|1|) does not 
necessarily hold in RPWIA. As we shall see in detail, this difference comes from the fact 
that in the non-relativistic approximation the nuclear current is expanded on a basis of 
free (positive energy) Dirac spinors, while in the relativistic approximation the nuclear 
current is directly written in terms of the relativistic bound nucleon wave function that 
contains both positive and negative energy projections in the complete Dirac basis. Thus 
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the question of whether factorization (in the sense of eq. ([]])) does or does not hold in 
RPWIA is intimately connected with the role of the negative energy projection of the 
bound nucleon wave function. 

Although the physical meaning of contributions from negative energy projections can 



be (and indeed is) a matter of debate [15], and their presence may be questioned, we would 
like to stress that in this respect our aim here is to analyze the quantitative importance 
of those contributions rather than to take sides on a particular conceptual interpretation 
or criterion. 

As we shall see the extent to which negative energy components are important depends 
on the choice of the current operator. Thus the present analyses is also useful in choosing 
a particular form of the current operator once one adhers to a given criterion. For instance 
if one would like to follow the point of view that negative energy projections should not 
contribute to physical observables one would choose a current operator that minimizes 
the role of such components. 

The organization of the paper is as follows. In section 2 we discuss in detail the relation- 
ship between RPWIA and PWIA differential cross sections. A separation of the RPWIA 
cross section and response functions is made into positive and negative energy projections. 
In section 3 we discuss various choices of the nuclear current operator and their influence 
on the positive and negative energy projections of the single-nucleon responses and cross 
sections. In section 4 we analyze the role of the negative energy projections and the lack 
of factorization in the total nuclear responses and differential cross sections. The main 
conclusions are summarized in section 5. 



2. Formal comparison of RPWIA to PWIA 

2.1 Differential cross section in RPWIA 

The general formalism for exclusive electron scattering reactions has been presented in 
detail in several previous papers. We refer in particular to refs. [Q and || for the nonrel- 
ativistic and relativistic treatments, respectively. Here we just summarize the kinematics 
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and focus on those aspects of the cross section that are of relevance to the points of 
discussion in this paper. 

In the (e,e'p) process an electron with four momentum = (e, k) is scattered through 
an angle 8 e to four momentum K' 11 = (V, k'). We denote the hadronic variables by P% = 
(Ma, O), Pg = (E B ,p B ), P£ = (E N ,p N ) the four momenta of the target, residual nucleus 



and outgoing proton, respectively. The target rest mass is Ma, and M = J E^ — p N 



Mb = y E B — p\ are the outgoing proton and residual nucleus masses (the latter may 
include possible internal excitation). As usual in these processes electrons are treated 
in the extreme relativistic limit (ERL): e = k, e' = k! . The four momentum transfer is 
given by = — K'^ = (uj, q). For the applications discussed in the next sections we 
consider 16 and have selected three kinematic situations: 

• I) q = 500 MeV/c, uj = 131.56 MeV 

• II) q = 1 GeV/c, uj = 432.8 MeV 

• III) q = 1 GeV/c, uj = 300 MeV. 

The value of uj in kinematics I and II corresponds to the quasielastic peak value (ujqe) 
while in kinematics III, uj < ujqe, the so-called y-scaling region. Denoting by p the three 
momentum of the struck nucleon (p = p m , the missing momentum, in the RPWIA, as 
discussed below), the range of variation of p-values has been chosen < p < 500 MeV/c. 
From similar arguments to those presented in ref. []TB[ it can be shown that the variation 



of the momentum is negligible, since one has that M B ^> p 2 . In fact, for A = 16 this 
variation is of the order of ~ 4% (kinematics I) and ~ 1.5% (kinematics II, III) in the 
whole range of p-values considered. 

In order to compare the content of the relativistic and non-relativistic treatments in 
plane-wave impulse approximation we focus, without lack of generality, on the differential 
cross section leading to a specific final state of the residual nucleus, corresponding to 
proton knock-out from a specific bound orbital b in the target nucleus. Following ref. [8], 
and taking unity spectroscopic factors and plane waves for incoming electron and for 
outgoing electron and proton, we write the differential cross section in RPWIA as 



6 



*o 2o? (s>\ PN MM B ^ M? (4) 



dQ e de'dn N g 4 \e J M A f r 
where X) denotes average over initial and sum over final polarizations. The transition 
amplitude is the contraction of the electron and nucleon («/$•) currents 

m = (5) 

with 

f e = K f (k')YKS k ) (6) 

J& = u„ N (PN)j&K b (p) , (7) 

where is the nucleon current operator to be discussed in section 3 and the momentum 
and energy of the bound nucleon satisfy, p = p N — q = —p B ; ej = — Ma + Mb + M = 
uj — T/v — Tb- The term denotes the Fourier transform of the relativistic bound 

nucleon wave function 

K b iP) = / dpe-V r K b (r) (8) 

with ty™ b (r) a solution of the Dirac equation with S-V potentials. For details see appendix 



A and refs. El 131. 



2.2 Positive and negative energy projections of the bound nucleon 
wave function and nucleon current 

It is known that only free solutions can be expanded in terms of free positive-energy 
Dirac spinors u alone. For an interacting (i.e., not free) relativistic wave function there is 
always a coupling to the free negative-energy Dirac spinors v. This implies that the free 
relation between upper (u) and lower (d) components 

a ' P (9) 



E -\- M 

does not hold in general for a bound relativistic wave function. As we shall see in detail, 
this is what differentiates RPWIA from PWIA. To analyze this point we insert in eq. (0) 



the completeness relation [117 



K(f> s)up(p, s) - v a (p, s)v p (p, a)] = 5 a/3 (10) 

s 

and write 

J]V = (Jn)u ~ (Jn)v > (11) 

where the first term, indicated by the index u, comes from the positive-energy projector 
involving the Dirac spinors u(p, s) and the second term, indicated by the index v comes 
from the negative-energy projector involving the Dirac spinors v(p,s). These contribu- 
tions are given by 

(r N } u = (J% = ^u(p N ,s N )>u(p,s) [u(p,s)^(p)\ (12) 

s 

(J&) v = (J») v = J2u(p N ,s N )J»v(p,s) Wp,8)*p(p)] . (13) 

s 

To simplify the notation, here and in what follows we suppress the index iV on nucleon 
currents and current operators. 

Using the explicit expression of the relativistic bound nucleon wave function in mo- 
mentum space given in appendix A, the u and v contractions in equations and (p~3| ) 
are found to be 



[u(p,s)^r(p)] = HY\j^^a Kb (p)(s\^(p)) (14) 



[v(p,s)K»(p)} = (-iY]j^^f3 K M(^-k(p)} > ( 15 ) 

where (s|$±£ fc (p)) indicates spin projections of the bispinors &± b Kb on a spin state ||s), 
and the radial functions in momentum space a Kb and f3 Kb are given by 

a Kb (p) = g Kb (p) - -g M S K J Kb {p) (16) 
A*(p) = ^^_ M 9K b (p)-S K J Kb {p) (17) 

with g Kb and f Kb the Bessel transforms of the standard upper and lower radial functions 
of the bound nucleon wave function in coordinate space (see appendix A for details). 

Obviously for a free nucleon /3 K (p) = 0. For a bound nucleon fl Kb (p) ^ in general, 
unless we impose the relation in eq. (|9]) which amounts to imposing (3 Kb {p) = 0. 
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2.3 Tensorial decomposition 

In order to connect with the standard PWIA we perform the usual decomposition of the 
transition probability in eq. (0j) into a leptonic and a hadronic tensor 

2J2\M\ 2 = Vllu W^ , (18) 

where 77^ is the leptonic tensor 

V = E &* (19) 



i f 



whose explicit expression can be found for instance in refs. (I], [18j , and W^ u is the hadronic 
tensor 

w» v = — — Yl JilJV * ■ ( 2 °) 

2jb + 1 s N m b 

We recall that in standard PWIA the hadronic tensor factorizes into a tensor for elastic 
scattering on a free proton W Ml/ and the momentum distribution of the (non-relativistic) 
bound orbital N b (p) 

(wn PW iA = w ,MU (p,q)N b (p) (21) 

with 

q) = J2 [u(Pn, s N )J»u(p, s)]* [u(p N , s N )J v u(p, s)] (22) 

S,S N 

and 

N b (p) = 7 -^—Y,\€ b (p)\ 2 , (23) 

where 0™ 6 (p) is the (non-relativistic) wave function of the bound nucleon in momentum 
space (see for instance refs. [[I], [Tj], [T9|j). 

The connection between RPWIA and PWIA can be best seen by substitution of 
eqs. ( |Il"l - |T5|) into eq. fl20|) . This leads to the following expression for the hadronic ten- 
sor 

= Wp v + + W£ u , (24) 

where Wp U (Wjff') is the contribution from positive (negative) energy projections only, 
while Wq 11 is a crossed term containing products of both positive and negative energy 
projections. These components of the hadronic tensor are given explicitly as follows 

^~EE(>):^)« (25) 

Z Jb 1 1 m b s N 

9 



wr 



N 



^Jb ' L m b s N 



= — -ee ((jr u (J u )v + (jrM u )u 

Using the relations 

2j b + 1 



m b 



5 XX * - 



8vr 



(s'\(T ■ p\s) 



III,, 



2]b + 1 
8pn 



it is easy to carry out the sum over m b in eqs. (|25| ~ p7|) to get (see appendix B) 



= N uu (p)W^ 

W^ = N vv (p)Z^ 
W£" = N uv (p)N^ , 
where we use the fact that in our case a Kb {p) and /3 Kb (p) are real to write 

N uu (p) = (a Kb (p)) 2 

n vv (p) = (Up)) 2 

N uv (p) = -2a Kb (p)(3 Kb (p) 



with 



MP) 



'E + M 
8vrM 
l E + M 



a« 6 (p) 



v 8vrM ' Kb 

and the single-nucleon tensors W^, Z^, M^ u are given by 



W" = Tr 



+ M\ T f ? N + M \ JU 



2M 



2M 



Z 



[IV 



Tr 



Tr 



2M 



- t ( £n±m; 

2M 



2M 



p 2M 
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where we use the notation J M = 7ot/ M+ 7o- Explicit expressions for these tensors for 
different current operators are given in appendix C. 

As already mentioned is the usual single-nucleon tensor appearing in the standard 
PWIA for unpolarized scattering, while and M^ v are new single-nucleon tensors that 
appear only when the bound nucleon wave functions contain a non-zero negative-energy 
projection. Using the identity 

7 ■ p = a ■ p7o7 5 (41) 



and following the arguments given in ref. [20] in the context of PWIA with polarized beam 



and target, it can be also shown that the single-nucleon tensor M^ v can be related to a 
diagonal tensor constructed from spinors quantized with respect to a spin axis pointing 
along a generic direction, TZ^ u (9r,(Pr), as 

N^ v = cos 0^"" (0,0) + sin# (cos 07^(^,0) + sin07^(^, ^)) (42) 

with 6, <fr defining the direction of the bound nucleon momentum p and 

1 



L J"(r N + M)J" (43) 



which is linear in the bound nucleon spin four-vector Sy\ The angles 9r, (f)R define the 
direction of the spin in the frame in which the bound nucleon is at rest. The explicit 
expression of 1Z^ U can be found in appendix C. Obviously this is not the only possible 



way to compute eq. ([!(]) — we have used eqs. (|2|, ^3|) because explicit expressions for the 



TZ^'s with the current operators that we consider here were already available. 

2.4 Comparison of RPWIA and PWIA 

In PWIA one defines a free electron proton cross section a ep 

2a 2 e' 



a 



ep 



Q 4 e 

to write the differential cross section as 

/ d 5 a \ p N MM B 



Vp»W^ (44) 



dn e de'dn N J PWIA M A f rec aePNb{p) (45) 
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with Nb(p) the momentum distribution of the non-relativistic bound orbital (eq. (f23|)) 
normalized to 1 (/ dpNb(p) — 1). 

In analogy, in RPWIA we may define single-nucleon cross sections o ep u , o~ ep v and o ep v 
corresponding to the single-nucleon tensors , Z^ v and appearing in the Wp U , 



W^ u and Wq" hadronic tensors: 



a ep 

uu 



a ep 



a ep 

" uv 



2a 2 


e' 


Q 4 




2a 2 


e 1 


Q 4 


£ 


2a 2 


e' 


Q 4 e 



r]^ v = a ep (46) 

(47) 

■ (48) 



Using the above definitions together with eqs. (|l"8D, (f24j) and (|30|)-(|32|) we can write 
the differential cross section in RPWIA (see eq. (§)) as 

= KlN ' M + aM) + " ZN ' M] • (49) 

where a^ u is the free electron-proton cross section appearing in eqs. (|1|) and ( f45|) obtained 
from positive-energy projections, while a^ p and cx^° are new components that are solely 
due to the negative-energy projections and that may only appear in scattering from a 
bound nucleon. Thus, in RPWIA the differential cross section depends on both positive 
and negative energy projections of the relativistic bound nucleon wave function. Equa- 
tion (f49|) is self-explanatory: it shows that one can perform a decomposition into a part 
proportional to the square of the positive-energy projection, N uu , with proportionality 
factor a ep u . This proportionality factor involves only free w-Dirac spinors and it is identical 
to the standard o ep that appears in the non-relativistic treatment (see eq. ([!])). The rest 
of the cross section is proportional (quadratically and linearly) to the negative-energy 
projection (N vv or N uv ) and the proportionality factor involves (quadratically or linearly) 
free v-Dirac spinors. 

From eq. ([H|) it is also a simple matter to go back to the non-relativistic limit by 
imposing condition @ on the bound nucleon wave function. Indeed if we impose 
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we get 



k = P { ^=0 (52) 
and therefore all the terms containing negative-energy projections in eq. (HO) become zero 



, 0) 2M g K 
a K = a y "> = \ = == (51) 



N<$ = (53) 
N<$ = , (54) 

while the term depending only on positive-energy projections becomes 

(55) 

Thus, we recover eq. (|49"D with the non-relativistic momentum distribution Nb(p) replaced 
by N^J(p) in eq. fl55|). Actually to go to the non-relativistic limit we should neglect terms 
of the order (1 — -||) and therefore we write 

KuiP) = , (56) 

where K is a normalization factor taking into account that g K (p) is not normalized to 1 
while N™u'{p) should be normalized to 1. Hence 

K- 2 = [ p 2 dpgl(p) . (57) 



We note that K > 1 because of the normalization condition of the relativistic bound 
nucleon wave function (see appendix A). Therefore in the non-relativistic limit the mo- 
mentum distribution N££'(p) > g%(p) / Air . Actually since f 2 (p) is small and does not con- 
tribute by more than a 3% to the normalization (for the case considered here K = 1.025) 
one has that (p) w Nj® (p) in the whole momentum range considered. 

In Figure 1 we show the functions N uu , N uv and N vv (in fm 3 ) corresponding to the shell 
lpi/2 in 16 0. The bound state wave function for the proton has been computed within 
the framework of the Walecka relativistic model. The mean field in the Dirac equation is 
determined through a Hartree procedure from a phenomenological relativistic Lagrangian 



with scalar and vector (S-V) terms. We use the parameters of ref. ||23[| , and the TIMORA 
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code [[12] to get the radial functions g Kb , f Kb in coordinate space and then transform to 
momentum space according to the equations in appendix A. 

In the left panel of fig. 1 we show the results for the three projection components of 
the momentum distribution N uu {p) (solid), N uv {p) (dotted) and N vv (p) (dashed) as given 
by eqs. (33 35). Note that the components N uu and N vv are positive whereas N uv can 
be positive or negative depending on the p-value. This is specified by the signs (+)/ (— ). 
The component N uu (p) clearly dominates in the region p < 300 MeV/c. In fact, in the 
maximum (p ~ 100 MeV) it is one order of magnitude larger than the N uv component 
and more than two orders of magnitude larger than the N vv component. This suggests 
that in this region one can expect negligible contributions from the components N uv (p) 
and N vv (p). Therefore, for low p-values, one may expect that projecting the bound 
nucleon wave function over positive energy states gives basically the same result as a fully 
relativistic calculation. This assesment of course needs some care because in order to get 
final results one also needs to evaluate the single-nucleon components <r^, a^ v and <r^. 

In the region of high missing momentum, p > 300 MeV/c, the situation is clearly 
different. In this p-region the values of the components N uv (p) and N vv (p) are similar to 
or even larger than that of N uu (p). Therefore, for p > 300 MeV/c one may expect that the 
effects of the dynamical enhancement of the lower components in the bound relativistic 
wave function will be observable. 

It is also interesting to compare the positive-energy projection of the full relativis- 
tic bound wave function with the non-relativistic limits defined in eqs. ( f^|56l) . In the 
right panel of fig. 1 we compare the positive-energy projection N uu (p) (solid) with N^(p) 
(eq. |55|) (dotted line) and with N™£'(p) ( ec L- HI) (dashed line). Note that the difference be- 
tween (p) and N™£- (p) is negligible for all p-values. Moreover, the difference between 
these two functions and the component N uu (p) is only visible for p > 400 MeV/c. Since 
up to p ~ 400 MeV there are no important differences between N uu and N™£', we shall 
speak roughly of positive-energy projected and non-relativistic limits when discussing 
total nuclear responses and differential cross sections in section 4. 
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2.5 Longitudinal- Transverse separation 

To finish with this section we note that one may also use the decomposition of currents 
into longitudinal and tranverse components and then write the different projections of the 
cr's as follows 

oZ = c7 MO tt{J2v k TZ^ (58) 
o% = ctmott (Y,v K K«j (59) 

°Z = °MOTT (J^VkKv) (60) 



with K = L,T, TL, TT the longitudinal, transverse, TL interference and extra tranverse 
contributions of the hadronic tensor for free nucleons. The single-nucleon responses TZ^ U , 
TZ^ V and 7Z^ V are given by taking the appropiate components of the single-nucleon ten- 
sors W M ^, Z^ u and J\f^ u , as is done for the R^ u in PWIA. It should be pointed out that 
contrary to the situation that occurs for on-shell nucleons, different results may be ob- 
tained depending on the choice of the current operator. In particular they depend on 
whether current conservation is or not imposed to eliminate one of the components, and 
on which of the 0- or 3- components is eliminated. This subject was discussed in detail 
in refs. p(| |21], |22| for the case of the TZ^ U responses with and without polarizations and 
will be discussed in the next section for the 1Z^ V and TZ^ V as well. Obviously, the hadronic 
response functions can be also written in RPWIA in the form 

R K = Rp + + Rc (61) 

with each component given by 

R K P = K« u N uu {p) (62) 

R% = K v N vv (p) (63) 

i£ = 72£X(p) • (64) 

In standard PWIA only the responses TV^ U which are due solely to the positive-energy 
components occur. In this case, by dividing the experimental 5-differential cross section by 
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cx^° and the kinematical factors, one could measure directly the momentum distribution 
N uu (p). However, in the more general case of RPWIA, as we see in eqs. ( f4T)| , |6"l~1) , this 
is not possible in general, as we have additional contributions from the negative-energy 
components that, a priori, cannot be isolated by experimental procedures. One could, 
however, study under what circumstances the results of eqs. (^,^l|) are dominated by 
some of the uu, vv or uv contributions. We present in detail this study in section 4 where 
results for the various observables of the scattering reaction are shown. 



Analogously to the analysis within the standard PWIA [£0], 21] the evaluation of 



the three projection components of the various single-nucleon quantities depends on what 
single-nucleon current operator is chosen, on whether current conservation is fulfilled or 
imposed, and in the latter case on how it is imposed. This study is presented in detail 
in section 3, where different prescriptions involving single-nucleon current operators and 
the continuity equation are considered. Here we simply mention that there exist some 
kinematics for which the uu, vv or uv contributions to the single-nucleon cross section 
and single-nucleon response functions are comparable, and therefore a fully relativistic 
analysis can give very different results from the standard PWIA. 



3. Current operators and single-nucleon responses 



3.1 Choices of the current operator for bound nucleons 



As discussed in several works |1|, |2(| ^T] |22| the choice of the current operator = J M is 
to some extent arbitrary. We discuss here the most popular choices denoted as CC1 and 



CC2 in ref. 21 



The current operator CC1 

J£ C1 = (F a + F 2 )Y -^{p + PnY (65) 
with P M = (E, p) (E = ^/p 2 + M' 2 and p = p N — q) is obtained from the current operator 



CC2 



J&n = FiY + ^^Qy ( 66 ) 
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replacing Q u by (P/v — P)v an d using the Gordon decomposition for the free nucleon case, 
i.e., assuming that initial and final nucleons in eq. satisfy the free Dirac equation 
(f — M)^ = 0. For free nucleons eqs. (|65|) and (^) are totally equivalent and, since 
the current is conserved, the time component can be eliminated by writing it in terms of 
the third component or vice versa. In the case discussed here the final nucleon in eq. (0) 
satisfies the free Dirac equation, but the initial bound nucleon does not, and therefore the 
nucleon current and the nuclear responses are different when using CC1 or CC2 operators. 
Moreover the current is not conserved and, once the current operator (CC1 or CC2) has 
been chosen, one may choose to impose current conservation by eliminating the third 
component (CCl^or CC2^), or the time component (CC1'- 3 *' or CC2^ 3 ^). One may also 
choose not to impose current conservation (NCC1 or NCC2). In general each of these 
six different choices gives a different result and it is not clear how to define criteria that 
favour one choice over the others (see for instance P^). In the absence of sharper criteria 
it is advisable to use as reference a current operator that guarantees current conservation 
for any initial and final nucleon wave function. For this purpose we choose the operator 

% = W + i^<rQ u -F 1 9jy. (67) 

for Q 2 < 0, which is also equivalent to the CC1 and CC2 operators in the free nucleon 
case. Although we define the operator only for the case of virtual photon (Q 2 < 0), 
we note that the last term in eq. (|6?j ) has no transverse component and therefore it does 
not contribute in the case of real photons anyway, where one could question its validity 
because of the divergence as Q 2 — > 0. This operator has several good features: a) with 
this operator the current in eq. (|7|) is conserved, and b) the amplitude M. in eq. (|5|) is 
invariant under the replacement of Jq by Jqc2 because the electron current is conserved, 
which implies that the last term in eq. fl6"7|) does not contribute to the amplitude M.. A 
similar operator was introduced by Gross and Riska in ref. p^ . 
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3.2 Comparison of different currents and single— nucleon responses 

To understand better the results shown in the next sections it is helpful to split the 
amplitude in eq. (|5p into its longitudinal (L) and transverse (T) components 

M=M L -M T (68) 

with 

M L = =j e °(j°-^/) (69) 

M T = ./;"./'• • #J» , (70) 

where we have taken the z-axis parallel to q, and used the continuity equation for the 
electron current {wj® = q ■ j e ). Now it is simple to check that the hadronic longitudinal 
current that we define for any current operator as 

Jt = (J° - ^ J3) (71) 

q 

is identical for the Jq and Jqc2 choices. Indeed 

J L - ( 7° - — 7 3 I - -— 7° - I 7° - — 7 3 I - J L (72) 
J c — \ J c q ) ~ q 2 \ CC2 J cc2 J — J cc2 ■ l' z J 

On the other hand since the transverse components are equal, 

Jcc2 = J c ' (73) 

both the longitudinal A4 L , and the transverse M. T amplitudes (and hence each of the 
response functions) are invariant separately under the replacement of Jc by Jcc2- We 
emphasize that this is so provided that the complete expression for Jq C2 in eq. is 
used. We note however that if one starts from the CC2 operator and imposes current 
conservation (as done in the CC2(°) and CC2^ 3 ^ versions), the longitudinal amplitude 
changes. In the CC2(°) version, one replaces J%c2 by ~Jcc2 an d therefore Jqc2 * s replaced 
by 



J&C2 ~^ Jcc2(o) — y 2 Jcc2 ■ (74) 
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Similarly in the CC2( 3 ) version one replaces Jqc2 by — Jcc2i or equivalently Jcc2 ls 

UJ 

replaced by 

JcC2 ~^ JcC2(3) = — ( ^2 J JcC2 • (75) 

In summary the NCC2 choice is equivalent to using the conserved current defined in 
eq. (|6?D . To the contrary, with the CC2^ or CC2^ 3 ^ choices the j\A. L amplitude becomes 
different and so do the response functions R L and R TL , while the transverse amplitude 
and response functions remain unchanged. 

The effect of going from the conserved current to the CC2^ or 002® choices can be 
easily evaluated by considering the differences 

^CC2(0) = ^002 ~ JcC2(0) (76) 



and 



^CC2(3) — JcC2 JcC2(3) ■ (77) 



As expected they are both proportional to ^> nQ^,^ with different proportionality factors, 



uj 



2 



Acc2 ( o) = ~^F^ N ^% (78) 

Acc2( 3 ) = Fi^n^ % ■ (79) 

Similar relationships hold when the CCl^ and CCl^ currents are compared to NCC1, 
but in this case the differences A(jci(o) anc ^ ^c*ci(3) a ^ so contain a term proportional to 
(1 — uj/uj). Indeed when comparing Jccito) an d JccifS) ^° ^cci one S e ^ s 



— UJ 2 

^CCl(O) = JcCl ~ -7(701(0) = ^ N ~^2 



(Fi + Fa) 2_^+^( 1 _? 

uj 2M V uj 



^C*C1(3) — JcCl "7(701(3) ~ 



uj 2M V uj 



(80) 
(81) 



Obviously the transverse currents satisfy 



J CC1(0) — J CC1(3) — J CC1 

jT _ jT _ rT 

J CC2(0) — J CC2(3) — J CC2 
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Finally we compare NCC2 to NCC1. We note that Jqc2 can be written without lack 
of generality as 

J&ci-J&ci = ^[{(M-r N )Y + r(M-f)} + (E + E N -E-E N )6 l 



2M 

A 

2M 



fj.,0 

{(M - f N )-f + r(M-f)}- 5^{u - aJ) 7 °f ] , i = 1, 2, 3 . 

(82) 



From eq. Q52] ) it is easy to see that all of the components are different for the NCC2 
and NCC1 choices when is not a free li-spinor 

J(JC2 Jccn ^ = 0,1,2,3 (83) 

and 

Jcc2 Jcci (84) 

Jcc2 7^ Jcci • (85) 

We shall discuss these differences in more detail in the next subsection. However it is 
also interesting to compare the positive-energy projections of J C C2 an d Jcci (^ ne w-parts 
defined in section 2). For the positive-energy projections one sees that 

(Jcc2)u = (Jcci) it (86) 

and 

Fn 

(J l cc2)uu = (J l cci)uu + 2^( U ~ uJ)un1°Yu , i = 1, 2, 3 . (87) 

Hence, when considering positive-energy projections one has that the time components are 
equal and the transverse and longitudinal components differ by an amount proportional 
to u — To. This in turn implies that 

(Jcc2(o))u = (Jcci(o))u j (88) 

while 

(Jcc2)u - (Jbci)u oc (oj - uJ)/2M (89) 

and 

(JcC2)u — {JcCl)u = (JcC2(0))u—(JcCl(0))u = (^CC2(3))« — (^CC1(3) )« x ( CJ— Cj)/2M . (90) 
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The above relationships are important in understanding how the different single- 
nucleon responses and the single-nucleon cross sections behave depending on the choice 
of the current operator and depending on whether we consider the positive-energy or 
negative-energy components. We show the behaviour of the various response functions 
for three different kinematical situations in figures 2-4. Figs. 2 and 3 are for kinemat- 
ics I and II, respectively (see section 2). As seen in figs. 2a to 4a the longitudinal (L) 
and transverse-longitudinal (TL) responses tend to change for each of the current choices 
(NCC1, NCC2, CC1(°\ CC2(°), CCl® or CC2( 3 >), while the transverse responses (T and 
TT) only depend on whether the CC1 or the CC2 current operator is used (figs. 2b to 4b). 
This follows from the fact that the transverse current is independent of whether current 
conservation is imposed or not. 

Let us first discuss the transverse responses whose behaviour is somewhat simpler 
to understand. We recall that NCC2 is equivalent to using the conserved current Jq 
and therefore we use it as a reference. As seen in fig. 2b when CC2 (thick line) is 
replaced by CC1 (thin line) the positive-energy component of the transverse responses 
K^ TT changes little while the negative-energy components and 7^ TT differ by 

an order of magnitude or more. The reason for this is easily understood by looking at 
eq. When considering the positive-energy component the first term in this equation 

becomes zero ((f — M)u = 0) and the difference between the two curves for TZ^ U is 
proportional to (lu — uj)/2M (see eqs. (f37|), fl90D); moreover in the interference response 
function TZ^u the (u — u)/2M dependence tends to cancel. On the contrary when one 
considers the negative-energy components V^J T and 7^ TT the first term in eq. (|82| ) is 
maximum ((f + M)v = 2M) and the difference between CC1 and CC2 is maximized. 
Although quantitative details depend on the kinematics, the same qualitative behaviour 
is seen in figs. 3b and 4b. It is also interesting to remark that with the CC2 choice 
the negative-energy components are of the same order as the positive-energy ones in the 
whole p-region considered, while with the CC1 current (because of (M — f) terms) the uv 
and v v components may become much larger than the uu ones. Hence one can expect a 
stronger dependence on the negative-energy projections f3 K when using the CC1 operator. 

These trends are also observed in the longitudinal and longitudinal-transverse re- 



21 



sponses (figs. 2a to 4a). The CC1 current gives similar responses as the CC2 current for 
the uu components, but for the uv and vv components, the CC1 current tends to give 
much larger results (for the absolute values). In addition, a very important role is played 
by the way in which current conservation is imposed (or not). If one chooses not to im- 
pose current conservation (NCC1 or NCC2) one observes the same qualitative behaviour 
as that just discussed for the transverse components. The choice of imposing current 
conservation by elimination of the third component (CC1'- *' and CC2^ )) also shows a 
behaviour similar to that discussed above and the results tend to be close to the NCC1, 
NCC2 responses. Particularly in this case CC1'- *' and CC2^ give identical results for 
the uu component (as is clear from eq. (|5^)). However the CCl^ and 002® choices 
cause large deviations in all of the components including the positive-energy components 
IZ^u 1 -, particularly for kinematics II. This behaviour is mainly due to the ^iQ^f term 
in eqs. (|T9|j8T| ) that (as seen in figs. 2a to 4a) may cause a strong reduction of the 



responses when one goes from NCC1 to CC1 (3) or from NCC2 to CC2^. This reduction 
is also present in the CCl^ and CC2^ choices, but in these cases it is largely attenuated 
by the uj 2 jq 2 factor in eqs. ( |78| , pOD and therefore the latter choices give results much closer 
to the NCC2 result that we use as reference. 

It is important to remark that the simple relations stated in eqs. (78 -[Sip are also 



responsible for the large differences observed between results obtained with CCl^ and 
CCl^ choices (or with CC2(°> and CC2( 3 )) in the context of PWIA (see refs. 0, H). 



In these references the study was restricted to a ep (or its various polarization components), 
although the cause of the large differences was not as clearly identified as it is here. In 
particular, the following interesting relationship between TL responses follows: 



-ryTL _ <t?TL 2 

<1, i = 1,2 . (91) 



"CCi /v CCi(0) _ U_ 
I^CCi ,K -CCi(3) 1 



-CCi (3) 

A more involved relationship can be established for the 1Z L responses obtained with dif- 
ferent choices. 

The single-nucleon cross section components cr^, <7^° and are also shown for 
completeness in figs. 5 and 6. The results shown correspond to a redefinition of the ex's 
including the factor M 2 /EE^ for easy comparison to previous work along similar lines, 
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i.e., 



M 2 



.ep 
uu 



ep 



(92) 



a. 



EE N 



with a ep as defined by de Forest [J2TJ. Similar definitions are used for a ep v and a ep v in terms 
of a ep v and <j^, respectively. 

Fig. 5 corresponds to kinematics I for 9 e = 30° (top panel) and 9 e = 150° (lower panel) 
and fig. 6 corresponds to kinematics II with 9 e = 12.5° (the forward-angle limit of the Hall 
A spectrometers at CEBAF). For the positive-energy components (er^) the dependence 
on the different choices of the current operator is maximal for kinematics I and forward 
angle. For kinematics II that dependence remains small at any angle. The dependence 
of o ep u (i.e., a ep ) on the current operator has been discussed at length in refs. pOj, |2~2fl . 
Much more notorious is the dependence of cx^ and <r^ on the various current choices. 
Obviously the behaviour of the single-nucleon cross sections at 30° in fig. 5 mainly reflects 
the behaviours of 1Z L (and partly of TZ TL ), while the behaviour of a's at 8 e = 150° (bottom 
in fig. 5) mainly results from the interplay between TZ TL and TZ T responses. Compared 
to the wild variation of a^ v and cx^, a^ u shows a rather mild dependence on the choice of 
current operator. 

We would like to emphasize that the changes in a ep u are minor for NCC1, NCC2, CC1*- '' 
and CC2(°) and are important for CC1'- 3 '* and CC2^ 3 \ particularly at small q (kinematics 
I) . As one increases q (q > 1 GeV) the latter choices may not be that poorly behaved. 

It is also important to remark here that the single-nucleon cross sections associated 
with the negative-energy projections [o ep v and <r^) tend to be larger in absolute values 
when the CC1 current operator is used. Therefore one may expect a larger effect of the 
negative-energy projections in the total differential cross sections when using Jqqi, and 
therefore a larger deviation from the non-relativistic PWIA result. 



We now study in detail the problem of different current operator choices for nucleons 
satisfying the Dirac equation with S-V potentials 



3.3 The effect of S-V potentials 



(f -M-V-fo + S)^ = . 



(93) 
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In principle both initial and final nucleons could be considered as solutions of this equation 
with the same potential; however, in the present work to simplify the problem somewhat 
in discussing the situation where the final nucleon is at relatively high energies (and 
so roughly quasifree) here we concentrate on the case (RPWIA) in which the outgoing 
nucleon is a solution of the free Dirac equation. Thus for the final wave function we have 



M)tyf = , 



(94) 



while for the initial 

(f - M)*i = (V l0 - S)% . (95) 

To understand the effect of the S-V potentials in the initial state we first write the 
wave function in coordinate space 



iEt 



X 

\<t>d{x)) 



with 



where p = —i'V and 



<t>d{x) = -~ 



<7 -p<p u (x) 



E(x)+M(x) ' 



(96) 



(97) 



E{x) = E - V(x) = E{x) 
M(x) = M - S(x) = M(x) . 



(98) 



For S and V constant (nuclear matter) the solutions are still plane waves satisfying 

E = ^k 2 + Ml (k<k F ) , E = E-V , M = M - S , (99) 
while for S and V x-dependent the solutions satisfy the standard equations 



x 



(x) 



with 



7/ 

V 



-(E 2 — M 2 )(f) u — irj'(x)cr ■ x<f) d 
■(E 2 - M 2 )<f) d + iri(x)(T -x(f) u 



x dx ^ ^ 



(100) 



1 d 

x dx 

24 



(S-V) 



(101) 



We are interested in working in momentum space and thus write (see also appendix 



A) 



Up) 



(27T) 



3/2 



(102) 



MP 



(2tt) 3 / 2 7 ^ av 7 (2tt) 3 / 2 

We now write (either in p or r-space)[j 



with 



dxe 



-ip-x (- i(T ■ Vju(g)) 



(103) 









* = 














L \e+mJ 





5^ 



(104) 



(105) 



and 



cr ■ p 



- " E + M-' (106) 
The first term in eq. ( |104| ) is proportional to the positive-energy projection of \& (it has 
zero negative-energy projection), while the second term is a correction proportional to the 
negative-energy projection of ^. We use this decomposition to analyze the dependence 
on the dynamical enhancement (<p' d ) of the bound nucleon wave function when different 
choices of the current operator are used. 

We start as before from CC2 and CC1 and particularize to the present initial and final 
nucleon wave functions satisfying eqs. (|95| ) and (p4]), respectively. Using eq. (|82]) together 
with eqs. Q93|) and Q95D we find that 

F 2 



TV _ TM 

J CC2 J CCl 



2M 



Y(V^-S) + (lu-uJ)5^ 



(107) 



with uj — uj = E — E,E = y/p 2 + M 2 and the energy of the bound nucleon wave function 

(-E = E N — oj = Ma — Mg, neglecting the recoil energy of the residual nucleus). Hence if 
*In this section p is used both for operator (when writing in r-space) and as c-number (when writing 
in momentum space) 
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we write down explicitly the fj, — 0, 1, 2, 3 components we find that 

F 2 



J CC2 J CC1 



r - T 6 

J CC2 J CC1 



u CC2 J CCl 



2M 
F 2 ^ 



V f [V-S'f + (u-w)]% 



2M 



(108) 
(109) 
(110) 



For the longitudinal components that enter in NCC2 and NCC1 one gets: 

F 2 



1 L — l h 

J CC2 J CC1 



2M 
F 2 ^ 



f 



2M 



/ 



V - S 7 ° + (a; - 57) - -7 3 (V^7° - S) 

Q 

V - S + u - ul}(fi u + ^<x • q{V + s)<p\ 



CO 



2 a ■ q(V - S)(fi u + [V + S + u -57] 



'III) 
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Hence there are in general three types of terms in the difference between the longitu- 
dinal CC2 and CCl currents: 

• One proportional to oj — 57 that is present even after projection into the positive- 
energy sector. We recall that only in CCl'- *' and CC2 ( ^ - ) choices the longitudinal 
currents are equal for the positive-energy projections (see eq. (|89|)). 

• One proportional to the SMALL nuclear (V — S) potential acting on the upper 
component (fi u . 

• One proportional to the large nuclear potential (S + V) acting on the down com- 
ponent (fid- 

The latter contribution tends to enhance the role of the negative-energy components. 
In other words any dynamical enhancement cfi' d of the lower component of the bound 
nucleon wave function may appear unphysically augmented in the nuclear response due 
to lack of current conservation. Similar considerations apply to the difference between 
the transverse CC2 and CCl currents, 



JX JX 

J CC2 J CC1 



F2_ [ -<t x (V + S)<fi d 

2M \ a x { y_ S)(j)u J 



(112) 
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4. Dependence of the response functions and reduced 
cross sections on the negative energy projections 

In this section we present results of response functions and differential cross sections for 
lPi/2 proton knock-out from 16 0. The aim here is to make a quantitative analysis of the 
relative importance of negative-energy projection contributions for various choices of the 
current operator. 

The total hadronic response functions in RPWIA are shown in fig. 7. We have chosen 
kinematics I and the forms of current operators discussed in section 3. As already men- 
tioned in that section the transverse responses (R T and R TT ) depend only on whether 
the current operator CC1 or CC2 is chosen, while the R L and R TL responses depend also 
on the current conservation prescription. We recall that NCC2 is equivalent to using an 
exactly conserved current. 

As expected the longitudinal response R L is practically identical for NCC2, CCl( ) 
and CC2(°) (also for NCC1, not shown), while it is substantially different for CC2^ 3 ). It 
is interesting to note that, unlike in PWIA, in RPWIA the longitudinal response R L for 
CCl^ 3 ) is closer to the CCl^ than to the 002® cases. The shift upwards of R L in the 
CCl^ case is due to the negative-energy projection because, as seen in fig. 2a, the 
TZ^ V single-nucleon response is much larger than TZ^ U , compensating for the smallness 
of the negative-energy projection of the wave function (N vv ) (see fig. 1). Clearly if we 
were to consider only positive-energy projections, similar results would be obtained with 
CC1 (3) and CC2 (3) (as is the case in PWIA) — the term 1Z% V is roughly an order of 
magnitude larger with CCl® than with 002® causing a larger longitudinal response for 
the relativistic bound nucleon. In the longitudinal response function R L the contributions 
of the negative-energy projections are maximized for CCl®, are minimized for NCC2 and 
are also negligible for CCl^ and CC2(°). This is illustrated in the top panel of fig. 8 where 
Rp, Rj^ and RJq are plotted separately for various current choices. Here we use P, N 
and C to denote positive, negative and cross terms, respectively. In the CCl^ case Rp 
is only 80% of the total R L , while for the other choices Rp is more than 95% of the total 
(see also figs. 10, 11). 
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As seen in fig. 7 the R TL response is the one where the effects of different prescriptions 
for the nucleon current are the largest. The maximum of this response can change by as 
much as a factor of 2 (or 3) when one goes from CCl( 3 ) to CC2( 3 ) and CCl^ (or to CC2(°>). 
Here one can see with the help of the lower panel of fig. 8 (see also figs. 10, 11) that the 
negative-energy projections of the bound nucleon play a very important role through the 
crossed term Rq L , while Rjf remains negligible. This follows from the fact that in the 
p region under study, TZ^ is of the same order as TZ^ and much larger than TZ^ (see 
fig. 3a). As in the case of R L also for R TL the negative-energy contributions are minimal 
with the CC2(°) choice and maximal with the CCl® choice. With this latter choice the 
negative-energy contribution Rq L is dominant, contributing by ~ 75% to the total R TL , 
with CCl^ and CC2^ 3 ) choices R^ L and Rj> L contributions are approximately equal and 
even with the 002^°^ choice the Rq L cross term gives a sizable (~ 30%) contribution. 
Since for all the choices the positive-energy projections in the maximum are similar, the 
net R TL responses in RPWIA vary strongly from one choice to the other. Since relativistic 
effects are so important in this response, even in the plane-wave limit discussed here, it 
is likely that experimental data may reveal such effects. 

For the transverse responses R T and R TT the differences seen in fig. 7 between CC1 
and CC2 results are solely due to the negative-energy projections (see also figs. 9,10,11). 
Again, these contributions (i?^' TT , i?^' TT ) are maximal for the CC1 choice. However, in 
the case of R T their total contribution at the maximum is less than 10% for CC1 and less 
than 3% for CC2. Relatively much larger are the contributions of the negative-energy 
projections to R TT , although this response is very small compared to R T and it is difficult 
to isolate R TT from differential cross section measurements. 

Focusing on the choices CC1 (0) and CC2(°> that are closer to NCC2, in all cases we 
show in figs. 10 and 11, respectively, the total R K responses and their P, N and C terms 
to better appreciate the relative importance of the negative-energy projections in each 
response function. Clearly for the large responses R L and R T the difference between the 
total relativistic responses and their positive-energy projection is negligible. This in turn 
implies that for these responses one cannot expect to see large effects in going from PWIA 
to RPWIA, since also N uu N™£\ Hence the factorization limit of PWIA is practically 
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preserved by RPWIA in the R L and R T responses except when the CCl^ 3 -* prescription 
is used for the current. On the contrary, the small response R TT and the longitudinal 
transverse response R TL are very sensitive to contributions from the negative-energy 
components, because the contributions from positive-energy projections are relatively 
small. For these responses the factorization limit of PWIA is badly broken even when the 
more conservative CC2^ ^ choice is used (see fig. 11). 

The breaking of this PWIA factorization limit in RPWIA has important consequences 
in the differential cross section. This is illustrated in figs. 12-14 for different kinematical 
situations. Fig. 12 is for kinematics I at two different electron scattering angles (left 
9 e = 30°, right 9 e = 150°), while figs. 13,14 are for kinematics II and III at forward 
angle (8 e = 12.5°). We compare the fully relativistic results given by the prescriptions 
CC1(°) (thin-solid), CC2(°) (thick-solid), CCl^ (thin-dash) and CC2( 3 ) (thick-dash) to 
the non-relativistic PWIA limit (dotted). The PWIA results in these figures have been 
calculated from eq. ( |45|) with N b (p) as given by N££'(p) in eq. (H) and with the CC1 (0) 
choice. It is important to stress here that these PWIA results practically coincide with 
the positive-energy projection contribution of the total RPWIA differential cross section 
in eq. fl49|) . Indeed the difference comes only from the difference between N uu (p) (eq. ([33])) 
and iV" n r '(p) (eq. (0)), which is negligible, as seen in the right-hand panel of fig. 1. Hence 
the differences between the dotted and thin-solid (CC1'- ' 1 ) lines in figs. 12-14 are due 
to the contributions from the negative-energy projections. Differences with thick-solid 
(CC2(°)), thick-dash (CC2< 3 >), and thin-dash (CC1®) are also partly due to the effect of 
the negative-energy contributions and partly due to the different current operator. 

Clearly seen in these figures is the fact that there are two important effects of the 
negative-energy projections of the bound nucleon wave function: A) A modification of 
the strength of the peak at low p (p < 300 MeV), where data are used to determine 
spectroscopic factors (particularly important for kinematics I and for the prescriptions 
CCl® and 002®); B) A substantial modification of the shape of the differential cross 
section, particularly in the high p region. This latter effect is apparent in all of the 
different kinematical situations considered. For instance in kinematics I one sees that 
even in the case 9 e = 150°, where the various prescriptions give practically the same 
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results in the low p region and the RPWIA practically coincides with PWIA, the high p 
region (300 < p < 500 MeV) shows important relativistic effects whose size depends on 
the choice of the current operator. 

We also point out that at high p the cross section is significatively higher for the CC1 
choices. The results shown in fig. 12 reflect the behaviour observed in fig. 7. A detailed 
quantitative analyses of the role of the negative-energy projections on the differential 
cross sections for different current choices at any p-value and various kinematics observed 
in figs. 12-14 can also be performed using the results shown in figs. 1, 5 and 6. 



5. Summary and final remarks 

We have studied the relationship between relativistic and non-relativistic treatments of 
the plane-wave impulse approximation to A(e,e'p)B reactions by inserting the complete- 
ness relation into the relativistic transition nuclear current. This allows one to separate 
the RPWIA differential cross section into contributions from positive-energy projections 
and from negative-energy projections of the bound nucleon wave function. This sepa- 
ration, exhibited in eq. ([49D , clearly demonstrates that the factorization limit of PWIA 
expressed by eq. (|1|) (see also eq. ((4"5])) breaks down in RPWIA due to the presence of 
the negative-energy projection that although small is non-zero for a relativistic bound 
nucleon wave function. Typically, N uv is one order of magnitude smaller than N uu and 
N vv is two orders of magnitude smaller for p < 300 MeV. This is so for the lpi/2 shell in 
16 considered here, and we have observed a similar trend for various shells in different 
nuclei . 



If we consider only the positive-energy projection of the relativistic bound nucleon 
wave function we recover the PWIA factorized expression that relates A(e, e'p)B scattering 
to free electron-proton scattering through the elementary differential cross-section a ep = 
a ep u . However if we also take into account the negative-energy projection of the relativistic 
bound nucleon wave function other single-nucleon components a^ v and o ep v appear that are 
not present in the case of electron scattering from a free nucleon (or antinucleon). It is in 
this sense that we speak of lack of factorization in RPWIA. Obviously, the extent to which 
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factorization breaks down is measured by the relative importance of the contributions from 
the negative-energy projections. 

We have studied the relative importance of these contributions with different choices 
of the current operators and of the kinematics. A quantitative study has been presented 
of — both single-nucleon and total — response functions and differential cross sections for 
different choices of the current operator and of the kinematics, focussing on the case of 
a lpi/2 bound nucleon in 16 0. The main outcome of this study can be summarized as 
follows. 

The role of the negative-energy component in the total differential cross section is 
to cause a reduction at low p (p < 300 MeV) and an increase at high p (p>300 MeV). 
This affects both the value of the spectroscopic factors and the shape of the extracted 
momentum distribution p(p) || ||, [TO]] . 

The reduction at low p is generally small (less than a 15% for the NCC1, NCC2, CC1 {0) 
and CC2 < -°' ) choices) and depends mainly on the kinematics. The largest reduction is found 
at small q values (q < 500 MeV) and at forward angles for the current choice CCl^ 3 -*. The 
sizeable difference at low p and q (forward angles) observed between the predictions of the 
CCl^ -* and CCl^ 3 -* choices can be traced back to the difference between the longitudinal 



currents in eqs. ([80,81), and it is also present in the positive-energy projections. Actually 
a sizeable difference at these kinematics is also observed between CCl^ and CCl® results 
(as well as between CC2^ and CC2^ 3 ^ results) in the single-nucleon component a uu . 

At low p, the positive-energy projections are dominant and the RPWIA results do 
not differ so much from the PWIA ones; however at p > 300 MeV the negative-energy 
projections play a more important role and the differences between current choices are 
enhanced. Indeed the single-nucleon components a uv and a vv in general depend much 
more on the choice of the current operator (CC1 or CC2) than does the positive-energy 
projection a uu and this is mainly seen in the total differential cross sections at high p. 

As a general rule, the NCC2 and CC2^ choices minimize the effect of the negative- 
energy components, while all choices corresponding to the CC1 operator tend to enhance 
the contributions of the negative-energy components. The analyses of the nuclear current 
components in coordinate space show that, compared to CC2, the CC1 choices give more 
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weight to the dynamical enhancement of the lower components that appear multiplied by 
the large S + V potential. 

The response functions R TL and R TT are particularly sensitive to the combined effect 
of the negative-energy component and the choice of the current operator, because the 
cross term R^ L (-R^ T ) can be larger than the positive-energy projection Rj, L (Rj, T ). As a 
net result one can find up to a factor of 3 difference between predictions of various current 
choices. While there is basically no difference between NCC2 and CC2(°) predictions, the 
choices NCC1, CCl^ and CC2^ 3 ^ predict an R TL response that is about 40% larger, and 
the CCl^ 3 -* choice predicts a three times larger R TL response. These large effects can be 
traced back to the difference between the CC1 and CC2 current operators in eq. 
Clearly the (f — M) term that has no effect in the positive-energy projection Rj> L has 
an important effect in the crossed terms Rq and in the negative-energy term R% L . 
Obviously the (f — M) term also causes differences between CC1 and CC2 predictions 
in the Rc and Rn terms of other response functions, but it is in the case of the TL 
and TT responses where the net effect is maximum because of the small values of the 
positive-energy projections. 

The effect in R TL is particularly interesting because this response can be more easily 
measured than R TT and data on R TL are already available. Whereas for most of the cur- 
rent prescriptions (NCC1, NCC2, CCl^ and CC2(°>) and at low p values (p<400 MeV) 
the relativistic and non-relativistic results for the longitudinal and transverse responses 
differ at most by a 15%, in the case of the longitudinal-transverse {R TL ) response the rel- 
ativistic prediction is at least a 50% larger than the non-relativistic one. It is important 
to point out that there is experimental evidence of the fact that a different spectroscopic 
factor is needed for the R TL response than for the R L and R T ones. This experimental 



evidence follows from a non-relativistic analyses of the data [[26|, [2~?fl . Although consid- 
ering other effects such as two-body meson-exchange currents can explain part of this 



difference |28| , it is clear that relativity plays a very important role in understanding this 
response, and it will be very interesting to see to what extent such discrepancies may be 
related to the relativistic effects found here p5 . 
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Appendix A 

• Free Dirac Spinors 

We follow the conventions of Bjorken and Drell JT7]]. In particular, the free 
Dirac spinors are 



Hp, s) 



v{p, s) 



where E denotes E = \JM 2 + p 2 . 

1. Normalization 

u(p, s)u(p, s) 

2. Completeness relation 



! E + M 
2M 



'E + M 
2M 



X 
cr p 

\e+m 

l ap 

E+M 



X 



X 



-v(p, s)v(p, s) 



(113) 
(114) 



115) 



Yl s ) u p(Pi s ) - v a(P, S)vp(p, s)] = 5 af: 



Relativistic Bound Nucleon Wave Function 



(116) 



The relativistic wave function for the bound nucleon in momentum space is given 
by 



Kip) 



( 



9k(p) 



(2tt) 3 / 2 



V 



cr p 



117) 
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with S K — k/\k\, j — \k\ — \ and the quantum number t given by the relation 

i k for k > 

(118) 



£ 



k for k > 
—k — 1 for k < 



The function $™(Q p ) is given by 

= m\jm) = £<^|jm>>7(W = T,X s (s\$7) (H9) 



and satisfies the relation 

The radial functions g K and / K in momentum space are obtained from the respective 
functions in coordinate space 

- / r 2 drg K (r)j e (pr) (121) 

7T JO 

A(p) = y-/ o r 2 drUr) 3l (pr) (122) 

with j^Qor) the Ricati-Bessel functions and £ = /c, ? = « — 1 for k > 0; £ — \k\ — 1, 
£ = \k\ for k < 0. The radial functions in coordinate space <? re (r) and f K (r) satisfy 
the Dirac equation 

df K k 1 



(i^ « - 1 

dr r 



f K -[E-M-U s -U v ]g K (123) 
<? K -[£-M + [/ 5 -[/ y ]/ K (124) 



with Us (Uy) the scalar (vector) potentials that describe the target nucleus [12,13], 
and the normalization is 

/ r 2 dr (gl(r) + f K (r)) = J p'dp (gl(p) + /») = 1 . (125) 

Hence the relativistic (vector) density in momentum space normalized to 1 is 

Nlip) = l 2 E VTK = { -^±lM . (126) 
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Appendix B 



In this appendix we present in detail the algebra needed in order to evaluate the positive 
and negative energy projection contributions to the hadronic tensor (see eq. (|24|)). Let 
us start with the contribution from the positive-energy projections. It is given by (for 
simplicity we suppress the index b on the bound proton wave function and quantum 
numbers) 



= ^ ££<>):<>>« 

47 ~r - 1 - m s N 

= 27Tt££[*(p>«) K(pW Hp,*') K(p)) 

J 771 ss' 

r a i * r A 

£ \u(p N , sjv) J M u(p, s)\ \u(p N , s N )J u u(p, s') 



X 



(127) 



S N 



Introducing the tensor W^/ s defined as 

WZ = Y,Hpn,sn)J»u(p,s) 



u (Pn, s N )J u u(p, s') 



:i28) 



and using the result obtained for the coefficient [u(p,s) ty™(p)] as given in eq. flUf) , we 
can write 

2 (E + M\ 



a 



,(p)I 2 £w^£( s |^>*( s '|^> • (129) 



2j + 1 V 2M 

Now we can make use of the general relation given by eq. fl2"8p to get the result shown in 
eq. ©. 

In the case of the contribution from negative-energy projections, we have 

2 



47 ~r 1 m sjv 

= 27TI ^ E s ) *') *^(p)] 

J Tn ss' 

fA ~i * r A 

m(Pjv, s N )J^v(p, s)\ [u(p N , s N )J u v(p, s') 

S N 



(130) 



Introducing the tensor given by 



s's 



£ [«(Pjv, s N )J ,1 v{p, s)\ \u(p N , s N )J u v(p, s 

s N 



:i3i: 
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and using the results given by eqs. (|15| , |120| ), we can write 

2 fE + M' 



wr 



N 



2j + 1 V 2M 
2 fE + M" 



2j + 1 V 2M 



n 

58' 



l^(p)i 2 E^E( s l $ -.)*( 5 'l^) 

ss' m 

l&(p)l a i£<2£ 
xJ(^-p)xlE(^n*W> • 



(132) 



Using the relation given by eq. (|28|) , we get 



E + M 



l^(p)l a iE^E[ 



,8tt/ \ M i ■■ jr 

\ / 1 88 

where the sum over the index 5 is simply 

E[x s + ^-p)x 5 



xj \?-p)x 



XJ'KV-PJX 



Xs'K* 7 -P)X 



(133) 



(134) 



and the final result reduces to the expression given in eq. fl3~Tp. Note that the single- 
nucleon tensor Z^ v = Y^s^ss can be easily calculated using trace techniques (eq. (|39|) .) 

Let us proceed now to the evaluation of the cross term containing products of both 
positive and negative energy projections. Proceeding as in the previous cases and using 



eq. fl29|) we can write 



-2 



1 \ fE + M 
8tt 



1 



/ F cot 



XAv-PlX 



where we have introduced a new single-nucleon tensor defined as follows, 



2 s' U s = E [ u (Pn, s N )J^u(p, s)]* [u(p N , s N )J u v(p, ft 



(135) 



(136) 



From general properties of the 7 matrices and making use of the trace techniques [17 
we can express the tensor X^ s as 



s ' s ~ 8M 2 



\- Tr { 75 (f + M){5 SS > + 1^s>s)J\?n + M)r} , (137) 



where the pseudovector (p^ s , (see ref. |20| for details) reduces to the four-spin of the bound 
nucleon Si in the diagonal case, s' — s. 
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It is clear from the above expression that the off-diagonal contributions are purely 
symmetric and therefore one may write 



-L-s's — &s's + Wss'-A , 



138) 



where A^ v is antisymmetric under /j, v and real, whereas S^, a is symmetric under 
H <-> v, real for diagonal terms and in general complex for off-diagonal terms. These 
properties combined with the relation 



xti* -p)x s 



X- S [<T-P)X~ 



allow one to write finally, 



-1\ (E + M 



8tt 



M 



a 



.W.(p)E«? 



with 



s ' s AM 



Tr 



(139) 



(140) 



141} 



It can be proved that J2 SS ' T^7 S xt'(^~^~)x s is j us t a trace and can be simply written as 



P 



Tr 



J' 



pff N + M\ o7 p 



2M 



p M 



2M' 



(142) 



with the tensor A/"^ as introduced in section 2.3 (see eq. ([TOD). 



Appendix C 



In this appendix we give the explicit expressions of the single-nucleon tensors W^, Z^ v 
and TZ^ for the two current operators CCl and CC2 in eqs. (B5,BS|). 

The following expressions for the various single-nucleon tensors are obtained: 

• with Jqqi current operator 



M 2 W^ = {F 1 + F 2 f [ P^P V N + P U P» + 



fJ,V 



F 2 (F 1 + F 2 )-F 2 2 



2 8M 2 



{p + p N y{p + p N y (143) 



37 



M 2 Z^ = (Fx + F 2 f ( P^P V N + P V P» - (P + Pjv) ' 



[XV 



Ft ^ 



8M 2 



Q\P + P N y{P + P N y + F 2 (F 1 + F 2 ) (P»P» N - P^P M ) (144) 



WRT = (F 1 + F 2 ) 2 (S£P^ + SZP»-P N -S L gn 



+ 



F 2 



AM 2 

El 
2 



p N ■ s l (p + p N y(p + p N y 



(Fx + F 2 ) ( S£(P + P N y + SUP + P N y) (145) 



with Jq C2 current operator 



M 2 W^ U 



f 2 ( p»p» n + Fp; + ^) +f 1 f 2 (q-q <r Q " QU + Q " Q " 



+ 



_F[_ 

AM 2 



p n ■ Q (p*Q v + P V) + P-Q (P N Q U + W) 



Q 2 (P»P U + P'P") - 2M 2 - 



g 
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+ 0"" 1 2M 2 g 2 



g 2 g 



2P W • gp • g 



(146) 



M 2 Z 



2 -s^f 



F 2 P»P V N + P^P£ 



(P + P 
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Fi f 2 [ W + W + Q^ + W _ g . ( p + Pn)9 ^ 



It 

AM 2 



p N ■ q (p*q v + p V) + p • g (W + p^g M ) 



- Q 2 (P»P U + P V N P 
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(147) 



WRT = F 2 (S£P» + S V L P& - P N ■ S L9 n 
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+ 
+ 



F ± F 2 



F? 



(S£<y + SZ<y-2Q-S L g 



AM 2 



Pn ■ Q {S1Q U + S V L Q») + S L -Q (P»Q V + P^) 
- Q 2 (P&SZ + P»S£) - P N ■ S L Q»Q» + <T (Q 2 P N ■ S L - 2P N ■ QS L ■ Q)} 

(148) 

with A ■ B = A^B*. 

Finally, we show in table I the components of the spin four-vector Si in the laboratory 
frame. 
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X cos <p 
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7+1 

i 1 + tSi) cos< ^ 
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T 7+1 
|^COS0 
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Table 1: Components of the spin four vector Sl(9 R ,(f) R ). The notation 7 = ||; x 
' E jj sin 9n — ff sin 9; x' = c °s # = ig- cos — 2/t and k = ^ has been introduced 
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Figure captions 

Figure 1: Left panel: projection components of the momentum distribution (in units of fm 3 ): 
N uu (p) (solid), N uv (p) (dotted) and N vv (p) (dashed). Right panel: N uu (p) (solid), 
Nj$(p) (dotted) and N™'(p) (dashed) (see text for details). 



Figure 2a: Projection components of the single-nucleon response functions, TZ" U , TZ^ V and TZ^ V . 

Kinematics I have been chosen and results for the pure longitudinal (top) and in- 
terference longitudinal-tranverse (bottom) responses are shown. Thick lines corre- 
spond to prescriptions that use the CC2 current operator and thin lines correspond 
to the CC1 operator. Prescriptions shown are: NCC2 (NCC1) (solid lines), CC2(°) 
(CC1 (0) ) (short-dash lines) and CC2® (CC1 (3) ) (long-dash lines). 

Figure 2b: Same as fig. 2a, except that now for the pure transverse (top) and tranverse- 
transverse interference (bottom) single-nucleon responses. Note that here all CC2 
(CC1) prescriptions collapse into a single thick (thin) solid curve. 

Figure 3a: Same as fig. 2a, except that now for kinematics II (see text). 

Figure 3b: Same as fig. 2b, except that now for kinematics II (see text). 

Figure 4a: Same as fig. 2a, except that now for kinematics III (see text). 

Figure 4b: Same as fig. 2b, except that now for kinematics III (see text). 

Figure 5: Projections of the single-proton cross section as given by eqs. (f4*6|-|4"8|) . Results cor- 
respond to kinematics I and 0a? — 0° (in-plane). Forward- angle electron scattering 
(8 e = 30°) (top panel) and backward (6 e = 150°) (bottom) have been chosen here. 
The labels are as in fig. 2a. 

Figure 6: Same as fig. 5, except now for kinematics II and forward-angle electron scattering 

(0 e = 12.5°). 

Figure 7: Hadronic response functions for the lpi/2 shell in 16 and kinematics I. The labels of 
the various curves are as in fig. 2a. Note that for the two pure transverse responses, 
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R and R , all the CC2 (CC1) prescriptions collapse into a single solid (dashed) 
curve. 

Figure 8: Components of the nuclear responses, Rp, R^ and R^ (see eqs. (|61]-|64])). Re- 
sults are shown for longitudinal (top panel) and tranverse-longitudinal interference 
(bottom panel) response functions. The rest of labeling is as in fig. 2a. 

Figure 9: Same as fig. 8, except that now for the transverse (top) and tranverse-transverse 
interference (bottom) responses. 

Figure 10: Hadronic response functions for the lpi/2 shell in 16 and kinematics I. All of the re- 
sults correspond to the prescription CC1^°\ For each response we compare the fully 
relativistic result (solid) with the three components as introduced by eqs. (^-Q): 
Rp (dotted), i?£ (dashed) and R§ (long-dash). Some of the curves are multiplied 
by different scale factors. 

Figure 11: Same as fig. 10, except that now the prescription CC2(°) has been chosen. 

Figure 12: Differential cross section for the lpi/2 shell in 16 0. Kinematics I have been chosen and 
results are presented for forward and backward electron scattering angles: 9 e = 30° 
(left panel) and 9 e = 150° (right panel). We show the fully relativistic results given 
by the prescriptions (see text): CC2(°) (thick-solid), CC1 (0) (thin-solid), CC2( 3 ) 
(thick-dash) and CCl^ 3 ) (thin-dash). We also show the results corresponding to 
the non-relativistic PWIA limit (dotted line) (see text for details). 

Figure 13: Same as fig. 12, except that now for kinematics II and forward scattering angle 

(0e = 12.5°). 

Figure 14: Same as fig. 12, except that now for kinematics III. 
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Figure 7 
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Figure 10 
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